ELSEVIER 


Journal  of  Power  Sources  65  ( 1997)  159-171 


JIBIRAL OF 


25 

Performance  and  modelling  of  a  direct  methanol  solid  polymer 

electrolyte  fuel  cell 

K.  Scott a’*,  W.  Taama a,  J.  Cruickshank  b 

a  Department  of  Chemical  and  Process  Engineering ,  University  of  Newcastle,  Newcastle  upon  Tyne  NE1  7RU ,  UK 
b  Defence  Evaluation  and  Research  Agency,  Haslar,  Gosport,  Hants  P012  2AG,  UK 

Received  4  November  1996;  accepted  30  November  1996 


Abstract 

The  performance  and  modelling  of  a  direct  methanol  fuel  cell  based  on  a  solid  polymer  electrolyte  membrane  (SPE)  is  reported.  Two  sizes 
of  cell  are  used:  a  small  cell  with  an  area  of  9  cm2  and  a  large  single  cell  with  an  area  of  250  cm2.  The  fuel  cell  utilises  a  vapourised  methanol 
fuel  at  a  porous  carbon /Pt-Ru  catalyst  electrode.  The  performance  of  the  fuel  cell  is  affected  by  the  cross-over  of  methanol  from  the  anode 
to  the  cathode  through  the  polymer  membrane  and  this  behaviour  is  modelled.  To  evaluate  cell  performance,  mathematical  models  are 
constructed  which  describe  mass  transport  in  the  porous  electrode  structures  and  the  potential  and  concentration  distributions  in  the  electrode 
regions.  These  models  are  used  to  predict  the  cell  voltage  and  current  density  response  of  the  fuel  cell. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  uses  methanol,  in 
the  form  of  vapour  or  liquid,  to  generate  electrical  energy. 
Advancements  in  electrocatalysis  have  enabled  this  electro¬ 
chemical  combustion  to  carbon  dioxide  to  proceed  without 
the  formation  of  intermediate  products  [  1  ] : 
anode: 

CH3OH  +  H20  =  C02  +  6H+  +  6e"  ( 1 ) 

cathode: 

3/202  +  6H+  +6e~  =3H20  (2) 

The  thermodynamic  reversible  cell  potential  for  this  over¬ 
all  reaction  is  1.214  V,  which  is  very  close  to  that  of  1.23  V 
for  the  hydrogen  fuel  cell,  and  consequently  has  generated 
the  interest  in  the  DMFC  as  an  alternative  power  source.  Fuel 
cells  generally,  among  their  many  uses,  are  expected  to  fill 
an  important  role  in  the  replacement  of  the  internal  combus¬ 
tion  engine.  The  DMFC  has  several  advantages  which  suit  its 
application  to  transportation,  including  high  efficiency,  very 
low  emissions,  a  potentially  renewable  fuel  source  and  fast 
and  convenient  refuelling.  This  latter  feature  is  considered  to 
be  significant  in  the  planning  of  fuel  supply  systems  which 
could  be  based  on  the  infrastructure  already  in  existence  for 
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internal  combustion  engines.  Hydrogen-based  fuel  systems 
are  less  convenient,  requiring  much  more  sophisticated  stor¬ 
age  and  transportation.  Alternatively,  reformation  of  metha¬ 
nol  on  board  a  vehicle  is  possible  although  a  significant 
weight  and  volume  penalty  results. 

The  current  advantage  of  the  hydrogen  cell  is  that  hydrogen 
oxidation  at  the  anode  is  very  fast  and  consequently  the  per¬ 
formance  of  the  hydrogen  cell  is  better  than  that  of  the  meth¬ 
anol  cell.  The  major  loss  in  voltage  of  hydrogen  fuel  cells,  at 
low  to  moderate  current  densities,  is  due  to  oxygen  reduction. 
This  requires  four  electrons  in  comparison  to  the  two  elec¬ 
trons  required  in  hydrogen  oxidation.  For  methanol,  six  elec¬ 
trons  must  be  exchanged  for  complete  oxidation  and 
consequently  the  oxidation  kinetics  are  inherently  slower.  A 
factor  of  the  slower  kinetics  is  believed  to  be  ‘poisoning’  of 
the  anode  catalyst  by  intermediates  formed  during  methanol 
oxidation  [2].  Oxidation  of  the  poisoning  intermediates  to 
carbon  dioxide  requires  the  adsorption  of  an  oxygen  contain¬ 
ing  species  (e.g.  OH,  H20) .  Adsorption  of  these  species  does 
not  occur  substantially  until  potentials  well  above  the  open- 
circuit  values  are  used  [3].  Platinum  by  itself  is  not  suffi¬ 
ciently  active  to  be  a  commercial  methanol  oxidation 
electrocatalyst.  Consequently  significant  activity  has  been 
seen  in  the  promotion  of  methanol  oxidation  with  significant 
results  achieved  with  the  use  of  binary  catalysts,  notably  Pt- 
Ru,  where  the  second  metal  forms  a  surface  oxide  in  the 
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Porous  Carbon 


Fig.  1.  A  schematic  diagram  of  the  direct  methanol  fuel  cell  with  solid 
polymer  electrolyte. 

potential  range  for  methanol  oxidation  [  4  ] .  The  area  of  meth¬ 
anol  oxidation  catalysis  needs  further  performance  improve¬ 
ments  to  push  the  DMFC  towards  successful  commercial 
exploitation  and  research  is  ongoing  in  several  laboratories 
worldwide. 

The  DMFC  can  function  in  several  electrolytes;  alkaline, 
acid  and  proton  conducting  polymer.  Alkaline  electrolytes 
pose  the  problem  of  carbonation  and  research  efforts  have 
focused  on  acid  electrolyte  and  more  recently  solid  polymer 
electrolytes  (see  Fig.  1).  Much  of  the  research  on  SPE  sys¬ 
tems  has  used  Nafion®  membranes,  typically  117,  from 
DuPont.  The  direct  methanol  fuel  cell  based  upon  solid 
polymer  electrolyte  (SPE)  has  the  additional  advantage  of 
no  liquid  acidic  or  alkaline  electrolyte.  Recent  developments 
in  electrode  fabrication  techniques  and  better  cell  designs 
have  brought  dramatic  improvements  in  cell  performance  in 
small-scale  DMFCs.  Typically,  power  densities  higher  than 
0. 1 8  W  cm  "  2  are  achievable,  and  power  densities  higher  than 
0.3  W  cm“2  have  been  reported  [5] .  These  power  densities 
are  however  substantially  lower  than  those  obtained  with 
hydrogen  fuel  cells,  0.6-0.7  W  cm-2,  with  which  loading  of 
the  platinum  anode  catalyst  can  be  substantially  lower  (0.1 
mg  Pt  cm  2). 

To  date  an  essential  condition  for  the  successful  operation 
of  a  DMFC  is  the  use  of  a  pressurised  oxygen  or  air  supply 
to  the  cathode.  Another  important  factor  is  the  concentration 
of  methanol  in  the  water-methanol  mixture  fed  to  the  anode. 
At  concentrations  higher  than  around  2  molar,  the  cell  voltage 
declines  significantly  due  to  poisoning  of  the  cathode  electro¬ 
catalyst  by  methanol  that  has  permeated  through  the  SPE 
(Nafion®)  membrane,  i.e.  methanol  crossover.  We  have 
measured  permeation  rates  of  water  and  methanol  using  a 
membrane-electrode  assembly  (MEA)  and  membrane  alone, 
installed  in  a  small-scale  cell  assembly  under  typical  DMFC 
operating  conditions.  The  results  obtained  have  been  used  to 
explain  the  relatively  poor  performance  of  the  DMFC  when 
low  pressure  oxygen  and  when  high  concentrations  of  meth¬ 
anol  are  used  [  6  ] .  Thus  a  second  area  for  research  to  improve 
the  DMFC  performance  is  in  polymer  membrane  electrolytes. 
Recent  work  has  reported  the  use  of  polybenzimidazole  [7] 
and  perfluorinated  sulfonimides  [  8  ]  as  polymer  electrolyte 


membranes  as  a  means  of  reducing  the  impact  of  methanol 
crossover. 

The  research  reported  here  relates  to  a  programme  geared 
to  the  scale  up  of  the  DMFC  to  a  stack  with  a  power  output 
of  0.5  kW.  Data  presented  are  for  single  cell  performance  of 
the  DMFC  using  Nafion®  1 17  membranes  and  the  modelling 
of  the  cell  to  predict  the  cell  voltage  and  current  density 
characteristics. 


2.  Experimental 

Tests  on  the  DMFC  were  performed  with  two  cells,  one 
with  a  cross-sectional  area  of  9  cm2  and  the  second  with  a 
cross-sectional  area  of  250  cm2.  Each  cell  was  fitted  with  one 
membrane  electrode  assembly  (MEA)  of  the  appropriate 
size.  In  the  small  cell  the  MEA  was  sandwiched  between  two 
graphite  paper  gaskets  (see  Fig.  2),  with  serpentine  flow 
paths  cut  out  for  methanol  and  oxygen /air  flow,  and  then 
between  solid  graphite  blocks  with  fluid  entry  ports.  The  cell 
was  held  together  between  two  aluminum  backing  plates 
using  a  set  of  retaining  bolts  positioned  around  the  periphery 
of  the  cell.  In  the  larger  cell  the  MEA  was  sandwiched 
between  two  graphite  blocks  into  which  the  flow  paths  for 
methanol  and  air  were  cut.  The  flow  path  consisted  of  a  set 
of  small  rectangular  channels  ( 1  mm  by  1  mm)  arranged  in 
a  serpentine  flow  path  (Fig.  3) .  The  fuel  cells  were  used  in  a 
simple  flow  rig  (shown  schematically  in  Fig.  4)  which  con¬ 
sisted  of  a  Watson  Marlow  perilstatic  pump  to  supply  aqueous 
methanol  solution  from  a  reservoir,  and  a  Eurotherm  temper¬ 
ature  controller  to  heat  and  vaporise  the  methanol.  Air  was 
supplied  from  cylinders  at  ambient  temperature,  bubbled 
through  water  for  humidification,  and  the  pressure  regulated 
at  inlet  by  pressure  regulating  valves.  All  connections 
between  the  cells  and  equipment  were  with  PTFE  tubing, 
fittings  and  valves.  The  output  from  the  cathode  side  of  the 
cell  was  analysed  on-line  using  an  AI  Cambridge  GC94  gas 
chromatograph. 

MEAs  studied  in  this  work  were  made  in  the  following 
manner.  The  anode  consisted  of  a  carbon  cloth  support  (E- 
Tek,  type  ‘A’ )  upon  which  was  spread  a  thin  layer  of  unca¬ 
talysed  (ketjenblack  600)  carbon,  bound  with  10  wt.% 


Fig.  2.  Schematic  diagram  of  the  experimental  small  scale  DMFC. 
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Fig.  3.  Photograph  of  the  large  scale  DMFC  showing  the  serpentine  flow  beds. 


Fig.  4.  Schematic  diagram  of  the  experimental  DMFC  flow  circuit. 


Nation®  from  a  solution  of  5  wt.%  Nation®  dissolved  in  a 
mixture  of  water  and  lower  aliphatic  alcohols  (Aldrich). 

The  catalysed  layer,  consisting  of  50  wt.%  Pt-Ru  (2  mg 
cm”2metal  loading)  dispersed  on  carbon  (ketjen)  and  bound 
with  10  wt.%  Nation®,  was  spread  on  this  diffusion  backing 


layer.  The  cathode  was  constructed  similarly,  using  a  diffu¬ 
sion  layer  bound  with  15  wt.%  PTFE,  and  1  mg  cm”2  Pt 
black  (Aldrich)  with  10  wt.%  Nation®  as  the  catalyst  layer. 
The  purpose  of  the  uncatalysed  layers  was  primarily  to  pro¬ 
vide  a  flat  surface  for  the  catalyst.  The  electrodes  were  placed 
either  side  of  a  Nation®  117  membrane  (Aldrich),  which 
had  been  previously  boiled  for  1  h  in  5  vol.%  H202  and  1  h 
in  1  M  H2S04  before  washing  in  boiling  Millipore  water 
(  <  1 8  mfl)  for  2  h  with  regular  changes  of  water.  The  assem¬ 
bly  was  hot-pressed  at  100  kg  cm”2  for  3  min  at  135°C.  The 
resulting  MEA  was  installed  in  the  cell  after  pressing,  and 
hydrated  with  water  circulated  over  the  anode  at  96°C  for 
several  hours. 


3.  Mathematical  model  of  the  DMFC  MEA  assembly 

The  model  structure  of  the  DMFC,  shown  in  Fig.  5,  con¬ 
sists  of  a  flow  channel  cut  into  a  graphite  flow-bed,  through 
which  the  reactant  flows:  adjacent  to  the  channel  is  the  dif¬ 
fusion  region  of  the  electrode,  comprised  of  a  highly-porous 
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carbon  cloth  backing  layer  and  a  thin  layer  of  uncatalysed 
Nafion-bound  carbon  black.  Next  comes  the  layer  of  porous 
electrocatalyst,  followed  by  the  Nation  membrane.  A  similar 
structure  exists  on  the  other  side  of  the  membrane.  The  model 
for  the  DMFC  needs  to  account  for  changes  in  potential  and 
allow  for  the  transfer  of  methanol  from  the  anode  to  the 
cathode  and  its  effect  on  the  performance  of  the  cathode. 

The  model  for  diffusion  through  the  carbon  fibre  cloth 
adjacent  to  the  electrocatalyst  layers  is  presented  in  the 
Appendix  for  the  case  of  the  methanol  anode  side  of  the  cell. 
The  effect  of  diffusion  through  the  highly  porous  structure  is 
negligible  under  realistic  operating  conditions.  Similar  con¬ 
clusions  can  also  be  made  for  the  oxygen  cathode  side  of  the 
cell. 

The  major  assumptions  adopted  in  the  model  are  as 
follows. 

(i)  Due  to  the  high  thermal  conductivities  of  the  graphite 
and  aluminum  cell  components,  the  cell  temperature  is 
assumed  to  be  constant  and  uniform. 

(ii)  Transport  along  the  flow  channels  can  be  described 
by  plug  flow. 

(iii)  The  pressure  is  uniform  within  each  cell  compart¬ 
ment  —  any  pressure  drop  occurs  across  the  membrane. 

(iv)  Due  to  the  thinness  of  the  diffusion  region  of  the 
electrodes,  transport  in  this  region  is  not  considered. 

(v)  Due  to  the  high  electronic  conductivity  of  the  carbon 
substrate  and  graphite  flow-beds,  no  voltage  drop  is  consid¬ 
ered  to  occur  through  the  thickness  of  the  electrode,  or  along 
the  flow  channels. 

(vi)  Due  to  the  amount  of  water  in  the  anode  flow  chan¬ 
nels,  the  membrane  is  considered  to  be  fully  hydrated. 

(vii)  The  concentration  of  reactant  is  taken  as  the 
weighted  average  of  the  inlet  and  outlet  concentrations. 

(viii)  Electrode  kinetics  can  be  described  by  the  Tafel 
equation. 

(ix)  The  overpotential  caused  by  methanol  crossover  is 
directly  proportional  to  the  concentration  of  methanol  at  the 
cathode. 

The  object  of  the  model  is  to  calculate  the  overall  cell 
voltage,  which  can  be  written  as: 

^cell  £cell  ^7an  ^7cat  ^7xover  ^7ohmic  ( ^ ) 

where  £ceU  is  the  difference  between  the  half-cell  potentials 
of  the  anode  and  cathode,  at  the  reference  current  density  i0 , 
corrected  for  the  thermodynamic  effects  of  temperature  and 
pressure.  In  the  first  instance  we  consider  that  the  anode  and 
cathode  overpotentials  and  rjcat  are  described  by  Tafel 
kinetics  at  the  electrodes,  and  a  one-dimensional  potential 
and  concentration  distribution  is  calculated  within  the  thick¬ 
ness  of  the  catalyst  layers.  Ohmic  overpotential,  7?ohmic,  is 
calculated  for  the  resistance  of  the  membrane,  and  the  effect 
of  methanol  crossover,  i.e.  the  crossover  overpotential,  over, 
is  calculated  from  the  flux  of  methanol  through  the 
membrane.  The  appropriate  expressions  for  these  potential 
contributions  are: 


Vxoyev  ~  XJ MeOH 


(4) 

(5) 

(6) 


3. 1.  Methanol  permeation  through  the  membrane 


Permeation  of  water  and/or  methanol  through  a  Nafion® 
membrane  will  take  place  under  the  driving  forces  of  con¬ 
centration  and  pressure  gradients,  and  electro-osmosis.  If  we 
assume  Fickian  diffusion  and  a  linear  concentration  gradient 
through  the  membrane  of  thickness  tm  (i.e.  the  effective  dif- 
fusivity  is  independent  of  concentration),  then  we  can  write: 

D  c2Kk  A 

j=  --Ac— r-AP+—  /  (7) 

/m  tlF 

Assuming  that  the  permeate  is  entrained  in  the  carrier  gas 
flow  at  a  rate  proportional  to  c2,  we  can  write: 

j  =  kc2  (8) 


Combining  Eqs.  (7)  and  (8): 


’-esM 


A  \  /(  D  K  A  i 
+-/  /  1+— +— AP 
F )/  \  ktm  ktm  ) 


We  have  measured  permeation  rates  for  water,  methanol 
and  a  water-methanol  mixture  through  Nafion®  117  and 
calculated  the  values  of  k  and  K from  the  data  [  6  ] .  To  describe 
the  effect  of  methanol  flux  on  the  cathode  performance  we 
assume  a  methanol  coverage,  0 ,  proportional  to  the  methanol 
concentration: 


@&C2,  MeOH  (10) 

Assuming  that  the  methanol-covered  fraction  of  the  surface 
area  has  a  lower  free  energy  for  oxygen  reduction,  in  which 
case  the  overpotential  produced  by  methanol  crossover  is 
proportional  to  the  flux  then, 

^?xover  A7MeOH  (11) 

where  \ is  an  empirical  constant  to  be  determined. 

This  model  predicts,  then,  that  the  overpotential  due  to 
methanol  crossover  will  have  a  current- independent  term, 
affected  by  the  pressure  differential,  and  a  current-dependent 
term  (producing  an  iP-like  drop)  due  to  electro-osmosis  of 
methanol.  By  measuring  the  effect  of  pressure  differential  on 
flux  and  on  overpotential  (correcting  for  kinetic  effects)  we 
can  determine  a  value  for  x  and  then  estimate  AMeOH* 

From  the  results  of  the  correlation  of  methanol  permeation 
data  and  the  calculation  of  the  appropriate  parameters,  a  sim¬ 
ple  model,  based  on  linear  electrode  kinetics,  has  been  pro¬ 
duced  based  on  Eq.  (7)  which  describes  the  impact  of 
methanol  permeation  on  the  cell  current  density-voltage 
response.  The  model  data  correlation  and  experiment  are  in 
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good  agreement  [  6  ] ,  in  the  current  density  range  50-350  mA 
cm“2,  which  is  the  range  of  practical  operation.  At  low  cur- 
rent  densities  the  model  is  not  expected  to  agree  due  to  the 
assumption  of  linear  kinetics.  With  the  introduction  of  But- 
ler-Volmer  kinetics  the  model  gives  a  good  fit  over  the  cur¬ 
rent  density  range  up  to  350  mA  cm-2. 

3. 1. 1.  Porous  electrocatalyst  layer  model 

The  model  is  used  to  predict  the  current  distribution  within 
the  porous  electrode  caused  by  a  mixture  of  poor  mass  trans¬ 
port  (diffusion)  and  low  protonic  conductivity.  Each  elec¬ 
trode  region  is  described  by  the  same  model  structure  with 
different  parameters  used  in  the  appropriate  solution.  To  illus¬ 
trate  the  model,  an  oxygen  consuming  cathode  is  considered 
as  shown  in  Fig  5.  As  oxygen  diffuses  into  the  electrode,  it 
will  be  consumed  at  a  rate  governed  by  the  local  effective 
overpotential,  r\{z ),  and  the  local  concentration  c(z),  where 
z  is  the  distance  into  the  electrode.  Assuming  a  Tafel  rela¬ 
tionship,  the  local  current  density,  i(z) ,  can  be  written  as: 

i(z)  =  iac(z) / cQ  exp{r)(z)  / (3}  (12) 

where,  i0  is  the  exchange  current  density  at  the  reference 
potential,  c0  is  the  reactant  concentration  at  the  oxygen /elec¬ 
trode  interface,  and  (3  is  a  constant  related  to  the  Tafel  slope. 

Fig.  5  shows  a  representation  of  the  porous  electrode  of 
thickness  /,  with  the  oxygen /electrode  interface  at  z  =  0  and 
the  electrode/ membrane  interface  at  z  =  1.  Protonic  resistiv¬ 
ity  within  the  electrode  is  given  by  p,  the  effective  catalyst 
surface  area  per  unit  volume  is  y,  and  D  denotes  the  effective 
reactant  diffusivity  within  the  porous  electrode. 

If  we  consider  the  potential  drop  across  a  small  distance 
5z  of  the  electrode,  we  can  write  Ohm’s  law  as: 

z 

(t>(z  +  8z)-<j>(z)=p8zfyj(z)dz  (13) 

0 

where  cf>(  z)  is  the  potential  of  the  proton-conducting  phase 
within  the  electrode.  In  the  limit,  as  8z  ->  0: 


Now,  the  rate  at  which  the  flux  of  reactant  (gradient  of 
concentration)  is  changing  within  the  electrode  is  related 
directly  to  the  local  current  density,  and  can  be  written  as: 


(18) 


where  F  is  Faraday  ’  s  constant  and  n  is  the  number  of  electrons 
involved  in  the  reaction. 

Eqs.  (17)  and  (18)  fully  describe  the  variation  in  concen¬ 
tration  and  effective  overpotential  as  a  function  of  the  distance 
through  the  electrode.  The  boundary  values  relevant  to  this 
problem  are: 


c(0)=Co, 


dc 

d  z 


dr) 

dz 


~0,  r)(D  —Et 

o 


(19) 

(20) 


The  solution  of  these  Eqs.  can  be  simplified  by  adopting 
reduced  variables  as  follows: 


In  the  case  of  the  methanol-consuming  anode,  an  additional 
term  should,  in  principle,  be  considered  in  the  material  bal¬ 
ance  equation  describing  the  effect  of  methanol  crossover  by 
an  electro-osmotic  drag  term.  It  is  assumed  that  the  extent  of 
this  methanol  drag  is  defined  by  a  drag  factor  g,  the  moles 
of  methanol  transferred  per  mole  of  proton,  and  the  local 
methanol  drag  is  given  by: 
molar  methanol  flux, 


d$(z) 

dz 


(14) 


The  measured  overpotential,  £,  is  the  difference  between 
the  potential  of  the  electronic  conducting  phase,  <f)c  (assumed 
to  have  negligible  resistivity ) ,  and  the  membrane  </>(  /) ,  minus 
the  reference  potential  difference,  A  </>0 

£=&-<£(/)- (15) 

The  effective  overpotential  at  any  point  within  the  elec¬ 
trode  can  be  written  as: 


viz)  =  4>c-(l>(z)~A(f>0 


(16) 


Therefore,  by  differentiating  Eq.  ( 14)  w.r.t.  z,  we  obtain: 


d2y(z) 
d  z2 


=pyi(z) 


(17) 


=  ndtagi{z)  /  F  (24) 

The  change  in  the  local  methanol  drag  is  given  by: 
dNmeJdz  =  n«jFdi(z)/dz=(ndljF)i(z)y  (25) 

The  material  balance  then  becomes: 

^=^'(z)(1+n"dra8)  (26) 


Thus  in  the  general  solution  of  the  model  the  inclusion  of 
a  methanol  electro-osmotic  drag  term  modifies  the  term  M2 
defined  below. 

Overall  the  basic  equations  which  enable  the  calculation 
of  the  voltage  components  in  Eq.  (3)  are: 


=  MhTcF  exp(Atf3  rT)r) 


=  M2xcr  exp(M3,r77r) 
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subscript  r  refers  to  the  region  considered  ( anode  or  cathode 
electrocatalyst  layer) . 


—  ^  —  Vloc  —  ^loc,  reac 

z  =  -r  V  =  C  = - 

l  Et  c, 


o, reac 


M,-(  '**>  ) 

Er  r  W^PcanF) 


J  MeOH 


t, 


^m,MeOH 


m 


t 


m 


J_  _  Q 

Pr”  TMeOH  =  XMe0H^H20 

a. 


m 


3.1.2.  Method  of  solution 

The  method  of  solution  requires  that  the  models  of  the 
catalyst  layers  give  the  required  overpotentials.  In  this  a  value 
for  Tjcat  is  chosen,  and  I  is  calculated.  A  value  for  17^  is 
obtained  using  the  given  /.  The  values  of  E ^  rjohmic  and 
TJxover  are  then  calculated,  to  obtain  the  final  cell  voltage,  VccU. 

The  equations  describing  the  concentration  and  potential 
distribution  within  the  electrode  are  solved  numerically  using 
the  finite-difference  method  and  Newman’s  BAND  algorithm 
for  the  resulting  simultaneous  non-linear  equations  (using 
modified  NL3BAND.C  software  [9,10]).  To  obtain  a  pre¬ 
scribed  value  of  /,  iteration  using  a  Newton-type  algorithm  is 
employed: 


/'(£*)  = 


given 


I(Ek)  —E 
/'(£*) 
/(£*  + 5) -/(£*) 


where  8  was  typically  ~  1 0  4  X  E.  The  subscript  refers  to  the 
k A  iterate  from  a  guessed  value. 


3. 1.3.  Temperature  effect  on  cell  voltage 
From  thermodynamics: 

AS 
~  nF 

where  E  is  the  electric  potential.  Assuming  A  S  is  constant 
over  the  temperature  range  considered,  then  the  change  in 
cell  voltage  can  be  written  as: 


AT- 


AS 

nF 


For  the  methanol  oxidation  reaction: 
CH3OH + ^02  -»  C02  +  2H20 


we  can  calculate  from  standard  thermodynamic  data  (at 
25°C ) ,  where  the  methanol  and  water  are  liquid: 

A H=  -726.51  kJ  mol-1,  AG=  -702.36  kJmol-1, 

A5=  -0.081  kJmol"1  K"1 


and  for  the  gaseous  case: 

A H=  -709.33  kJ  mol-1,  A G=  - 889.42  kJ  mol" 1 , 
AS=  +  0.604  kJ  mol  - 1  K~ 1 
so  that: 

©  =  -0.140  mV  K-1,©  = +1.043  mV  K“‘ 

'  '  liq  '  '  gas 


3. 1.4.  Pressure  effect  on  cell  voltage 
From  thermodynamics: 

(dE\  _  AV 

[dpJT~  ~  nF 

which  gives,  on  integration: 


ANRT ^ 
AE= - In  — 


nF 


For  the  oxygen  reduction  reaction: 


2H++^02  +  2e 
2 


H,0 


we  get  AN  for  the  liquid  and  gaseous  cases  as  -0.5  and 
■F  0.5,  respectively.  At  80°C,  increasing  the  pressure  to  2  bar 
will  increase  the  voltage  by  10.5  mV. 


3.1.5.  Conductivity  of  Nafion 

This  can  be  described  by  the  following  relation  [11]: 


ref 


O' m  =  exp 


As  a  reference  point,  we  use  a.c.  measurements  performed 
in  our  laboratory  which  give  the  conductivity  as 
0.073  +  0.008  S  cm,  with  a  similar  value  for  recast  Nafion 
films,  at  25°C.  This  gives: 


crm  =  0.073  exp 

i 


S  cm 


3.1.6.  Diffusion  coefficients  of  water  and  methanol  in  Nafion 
117 

The  variation  of  diffusi  vity  with  temperature  can  generally 
be  described  by  an  equation  of  the  form: 


> 

1  lVI 

A  =  A,  ref  exp 

_  R 

KTrct  Tj  J 

For  the  diffusion  of  water  through  Nafion  1 155  E.W.,  Yeo 
and  Eisenberg  [12]  calculated  a  value  for  A E/R  of  2416  K. 
Springer  et  al.  [  13]  found  a  value  of  2436  K. 

We  will  use  the  latter  value,  as  it  pertains  to  measurements 
on  Nafion  117.  For  a  reference  value,  we  use  the  PGSE 
measurements  of  Zawodzinski  and  Springer  [14]  which 
give: 
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Dm  water  =  7.3  X  10  6  cm2  s  1  @80°C 
so  that: 


D 


m,  water 


—  7.3  X  10-6 


exp 


For  methanol  in  Nafion,  we  use  the  same  activation  energy 
as  for  water,  since  there  is  a  paucity  of  results  in  this  area. 
For  a  reference  point,  we  refer  to  the  work  of  Kauranen  and 
Skou  [  15]  who  measured  a  (superficial)  diffusivity: 

An.  MeOH  =  4-9X  10-6cm2  s'1  @60°C 


giving: 


D 


m.  MeOH 


=  4.9  X  10“6 


exp 


3. 1. 7.  Diffusion  coefficients  in  the  catalyst  layers 

The  diffusion  of  methanol  in  the  catalyst  layer  is  assumed 
to  have  a  similar  temperature  dependence.  Using  a  reference 
value  for  methanol  in  water  at  80°C  given  by  Kauranen  and 
Skou  [  1 5  ] ,  so  that: 

A,  MeOH  =  2.8  X  1(T 5  exp[2436^^  - 

When  multiplied  by  the  porosity-tortuosity  factor,  e,  we 
obtain  the  effective  diffusion  coefficient.The  value  of  AMeOH 
is  assumed  to  be  given  simply  by: 

AMeOH  "  XMeOH Ah20 

where  AH2Q  has  been  given  by  Zawodzinski  and  Springer 
[  14]  as  2. 0-2. 9  H20/H+  in  fully  hydrated  Nafion  1 17.  We 
will  use  the  average  value  of  2.5. 

3.1.8.  Dependence  ofi0  on  temperature 
Parthasarathy  et  al.  [16]  have  given  the  temperature 

dependence  of  the  electrode  kinetics  of  oxygen  reduction  as 

h,  cat  =  '0e  cat  exp 

For  methanol  oxidation  at  the  anode,  Troughton  [1]  has 
measured  the  activation  energy  for  a  Pt-Ru  supported  catalyst 
as  70  kJ  mol  -  \  giving: 

<0,  an  =  *0e  an  exp 

The  Tafel  slope  for  the  oxidation  reaction  can  also  be  esti¬ 
mated  from  Troughton’ s  work  as  46  mV  /  factor  of  e  at  80°C. 

3.1.9.  Concentrations 

A  methanol  concentration  of  one  molar  translates  into  the 
following  value  for  the  mole  fraction: 

Co,  MeOH  =  1  X  10-3  mol  cm-3 

*0,  MeOH  =  0 .0 1  84 


8804J— —  —  J 
\T„.f  T 


For  the  oxygen  concentration  at  the  cathode,  we  assume 
ideal  gas  properties,  so  that: 


co,  o2  ~ 


10  XRT 


mol  cm 


-3 


4.  Results  and  discussion 

One  purpose  of  the  model  is  to  estimate  the  variation  in 
local  reactant  concentration  and  potential  to  predict  the  total 
polarisation  potential  at  both  electrodes.  Figs.  6  and  7  show 
the  typical  distribution  of  potential,  concentrations  and  cur¬ 
rent  density  for  the  cathode  and  the  anode,  respectively.  The 
activity  of  the  electrode,  seen  in  terms  of  local  current  density 
and  overpotential,  is  higher  near  to  the  membrane  and 
decreases  towards  the  carbon  backing  layers  as  current  tends 
to  follow  the  more  conductive  electronic  path  through  the 
electrode  structure.  In  the  case  of  the  cathode,  at  what  is  a 
relatively  low  overall  current  density  of  approximately  62 
mA  cm-2,  there  is  an  approximate  25%  reduction  in  the 
oxygen  concentration  from  the  gas  side  towards  the 
membrane.  The  extent  of  the  overall  polarisation  of  the  elec¬ 
trode  depends  on  several  parameters  and  notably  the  thick¬ 
ness  of  the  electrode.  The  data  in  Fig.  6  are  for  a  relatively 
thick  oxygen  electrode  (200  pm) .  The  overall  potential  ver¬ 
sus  current  density  for  this  electrode  is  shown  in  Fig.  8.  The 
effective  Tafel  slope  for  the  electrode  is  approximately  twice 
that  for  a  non-porous  electrode,  i.e.  the  extent  of  polarisation 
of  the  electrode  is  increased  as  a  result  of  the  distribution  of 
current  in  the  porous  structure.  This  behaviour  corresponds 
with  that  predicted  by  simple  models  of  the  electrode  based 
on  an  assumption  of  a  constant  overpotential  in  the  electrode, 
derived  from  Eq.  (8).  This  predicts  for  thick  electrodes  (or 
high  current  densities,  poor  diffusion)  that  the  total  current 
is  given  by: 

I  =nFDC0yia  exp (  17/ 2/3) 

The  distribution  of  methanol  concentration  and  overpoten¬ 
tial  in  the  porous  anode,  as  shown  in  Fig  7,  is  for  a  high  value 
of  overall  current  density  at  a  thin  catalytic  electrode  (5  pm) . 
The  extent  of  polarisation  produced  by  the  porous  electrode 
is  relatively  small  even  though  there  is  a  large  change  in  the 
methanol  concentration  through  the  electrode  structure. 

The  complete  model  of  the  MEA  is  used  to  predict  the 
variation  of  cell  potential  with  overall  current  density,  both 
with  and  without  the  influence  of  methanol  crossover.  The 
experimental  work  has  shown  that  there  is  a  significant  effect 
of  increasing  the  oxygen  pressure  on  cell  performance  which 
cannot  be  predicted  from  thermodynamic  or  kinetic  behav¬ 
iour  or  by  the  model  which  does  not  allow  for  the  effect  of 
methanol  crossover,  as  shown  in  Fig.  9.  The  predictions  in 
Fig.  9  represent  the  ideal  fuel  cell  characteristics  if  methanol 
crossover  could  be  eliminated.  When  the  factor  of  methanol 
crossover  is  introduced  there  is  a  significant  change  in  the 
cell  characteristics  (see  Fig.  10),  cell  voltages  are  signifi- 
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Reduced  Distance  through  Electrode 

Fig.  6.  Variation  in  the  local  values  of  concentration  of  oxygen,  overpotential  and  current  density  in  the  oxygen  cathode.  M2  =  16,  M3  =  0.6,  Er  —  0A,  (1=0.06. 


Distance  through  Catalyst  Layer,  5  microns 

Fig.  7.  Variation  in  the  local  methanol  concentration  and  overpotential  in  the  anode.  /= 938  mA  cm-2.  Methanol  feed  concentration  -  1.0  kmol  m 


cantly  lower  and  decrease  when  the  oxygen  pressure  is 
reduced. 

The  agreement  between  the  model  and  the  experimental 
data  for  the  small  scale  cell  (see  Fig.  10)  is  good  over  a 
significant  range  of  practical  operating  current  densities,  up 
to  approximately  350-400  mA  cm'2.  The  poorer  agreement 
at  low  current  densities  is  partly  due  to  the  assumption  of 
Tafel  kinetics  and  the  relatively  simple  model  for  the  effect 
of  poisoning  of  the  cathode  due  to  methanol. 

The  experimental  data  shown  in  Fig.  10  are  for  a  relatively 
low  operating  temperature  of  80°C.  The  maximum  power 
densities  are  of  the  order  of  0.2  W  cm-2.  Power  densities 


reported  by  other  researchers  are  generally  at  higher  temper¬ 
atures,  for  example  Ren  et  al.  [  17]  report  values  of  0.38  W 
cm  “ 2  at  1 30°C,  using  oxygen  at  a  pressure  of  5  bar,  and  0.25 
W  cm-2  at  1 10°C  using  air  at  3  bar.  A  difference  of  50°C  in 
the  operating  temperature  makes  a  significant  difference  in 
the  cell  performance.  For  example,  at  a  constant  cell  voltage 
of  0.5,  Ren  et  al.  report  a  reduction  in  current  density  from 
around  450  mA  cm'2  to  about  170  mA  cm-2  when  the 
temperature  is  reduced  from  130  to  80°C.  In  addition,  these 
values  were  obtained  using  Nation  1 12  membrane  as  opposed 
to  the  Nation  117  used  in  this  work.  Nation  112  shows  approx¬ 
imately  half  the  cell  resistance  of  Nation  117  when  used  in  a 
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Current  density  /  mA  cm'2 

Fig.  8.  Predicted  variation  in  the  cathode  potential,  current  density  characteristics  for  the  oxygen  cathode.  M2  =  16,  M3  =  0.6,  £,.  =  0.1,  /0  — 0.01,  £>  =  0.1  cm2 
s_l,c0  =  3e'5  mol  m-3. 


Fig.  9.  Performance  of  the  DMFC  without  methanol  crossover.  80°C,  1 .0  M 
methanol,  oxygen  pressure  (gauge)  +0:0  bar,  +1:1  bar,  +  2:2  bar.  (  +0, 
+  1  and  +  2  in  Fig. ) . 

DMFC  [17].  When  these  two  factors  of  temperature  and 
membrane  material  are  allowed  for,  the  recently  reported 
performances  of  the  DMFC  are  very  comparable.  For  exam¬ 
ple  0.35  W  cm”2  at  97°C  [5]  and  0.385  W  cm"2  at  130°C 
[  17] ,  both  using  oxygen  at  5  bar. 

The  use  of  oxygen  in  the  DMFC  clearly  has  benefits  over 
air  in  terms  of  improved  power  performance.  However,  most 
applications  of  the  DMFC  will  require  the  use  of  air  as  oxidant 
at  a  modest  pressure.  This  factor,  as  well  as  issues  related  to 
safety  (in  scaling  up  to  an  eventual  0.5  kW  stack) ,  has  meant 
that  scale-up  studies  have  been  entirely  restricted  to  the  use 
of  air.  Fig.  1 1  shows  a  performance  of  the  DMFC  operating 
at  a  temperature  of  1 13°C  with  1.0  M  methanol  solution.  The 
maximum  power  density  is  1 84  mW  cm  "  2.  The  performance 
of  the  single  large-scale  cell  is  generally  comparable  to  that 
of  the  small-scale  cell  in  terms  of  power  output  and  I/V 
response.  However,  the  performance  is  not  identical. 


Typically,  we  see  a  somewhat  poorer  performance  with 
the  large-scale  cell  under  what  are  essentially  identical  meth¬ 
anol  feed  conditions.  Factors  that  can  explain  this  include  the 
reproducibility  of  preparation  of  MEA  for  the  larger  scale. 
Other  major  differences  in  behaviour  on  scale-up  are  changes 
in  the  fluid  temperatures  in  the  cell,  the  %  conversion  of 
methanol  achieved  and  factors  affecting  the  current  distri¬ 
bution  over  the  electrode  surface.  Investigation  of  these  fac¬ 
tors  is  the  subject  of  ongoing  research  into  the  scale-up  of  the 
DMFC. 

An  increase  in  scale  of  operation  of  the  DMFC  will  result 
in  changes  in  the  electrical,  fluid  mechanic  and  thermal 
response  of  the  system.  The  operation  of  a  cell  stack  contain- 


Fig.  10.  The  effect  of  methanol  crossover  on  the  DMFC  performance.  80°C, 
1.0  M  methanol,  oxygen  pressure  (gauge)  +0:0  bar,  +1:1  bar,  +  2:2  bar 
(  +0,  +1  and  +2  in  Fig  ).  Experimental  4:0  bar,  5:1  bar,  6:2  bar  (series 4, 
5  and  6  in  Fig.). 
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Fig.  11.  Performance  of  the  large  scale  DMFC.  1 13°C,  1.0  M  methanol. 


ing,  for  example,  25-50  membrane /electrode  assemblies 
capable  of  power  output  in  the  range  1-5  kW  will  require  an 
efficient  flow  and  thermal  control  system. 

Heat  generated  by  the  cell  stack  could  be  of  the  same  order 
of  magnitude  as  the  designed  electrical  power  output  and  will 
be  available  principally  in  the  exhaust  fuel  stream  and  also  in 
the  exhaust  air.  The  system  will  have  to  perform  the  following 
functions: 

•  supply  humidified,  pressurised,  preheated,  clean  air 

•  recover  water  vapour  and  methanol  (due  to  crossover) 
from  the  exhaust  air  stream 

•  supply  methanol  as  a  hot  aqueous  liquid  or  a  vaporised 
feed 

•  separate  carbon  dioxide  product  gas  from  the  exhaust  feed 

•  recover  the  heat  in  the  exhaust  feed  to  preheat  feed  and  air 
streams 

•  supply  methanol  to  the  cell  feed  to  maintain  the  desired 
optimum  concentration  for  cell  power  performance 

To  study  and  understand  the  interactions  of  the  above  func¬ 
tions  in  the  fuel  cell  system,  a  thermal  model  of  the  system 
itself  is  used  to  predict  the  temperature  changes  which  occur 
in  both  fuel  and  oxidant  streams  in  the  cell  as  a  function  of 
cell  operating  conditions  in  dynamic  operation  [18].  A  sec¬ 
ond  function  of  the  model  is  to  establish  the  level  of  internal 
temperature  changes  in  the  cell  which  have  implications  in 
the  mechanical  function  of  MEA  components.  The  compo¬ 
nents  of  this  model  consider  that: 

(i)  Catalyst  layers  are  thin,  compared  to  the  membrane, 
and  act  as  either  uniform  heat  source  or  sink. 

(ii)  The  bipolar  plates  have  an  overall  heat  transfer  coef¬ 
ficient  which  acts  between  the  two  separated  fluids. 

( iii )  The  membrane  has  a  high  thermal  resistance  and  Joule 
heat  is  liberated  into  both  fluids. 

(iv)  Heat  is  transferred  by  flow  (electro-osmotic)  of  water 
and  methanol  across  the  membrane  from  anode  to  cathode. 


( v)  Heat  transfer  through  the  porous  carbon  layers  is  taken 
into  account. 

This  model  of  the  MEA  will  give  the  required  heat  flux 
terms  to  use  in  dynamic  thermal  models  of  the  fluid  flow  in 
the  channels  of  the  flow  beds  incorporated  into  the  bipolar 
plates.  The  current  flow  beds  used  in  the  DMFC  are  made  as 
sets  of  small  rectangular  channels  (approx.  1  mm  by  1  mm) 
arranged  in  a  serpentine  flow  path.  This  therefore  means  that 
there  are  seven  possible  variations  in  the  overall  correspond¬ 
ing  flows  of  both  fluids;  co-,  counter-  and  cross-flow.  The 
above  thermal  model  is  solved  with  a  set  of  equations  describ¬ 
ing  the  change  in  pressures  and  material  balances  for  meth¬ 
anol  consumption,  oxygen  consumption,  water  production, 
and  membrane  transport  and  overall  mass  balances. 


5.  Conclusions 

The  results  of  this  work  have  shown  that  acceptable  per¬ 
formance  of  the  polymer  electrolyte  DMFC  can  be  achieved 
at  the  modest  temperature  of  80°C  using  vaporised  aqueous 
methanol  feeds.  The  importance  of  utilising  a  high  cathode 
oxygen  or  air  pressure  to  achieve  good  performance  is  also 
shown.  The  cell  potential  versus  current  density  characteris¬ 
tics  can  be  reasonably  well  predicted  using  a  mathematical 
model  of  the  MEA  based  on  the  variation  of  reactant  concen¬ 
trations  and  overpotentials  in  the  catalyst  layers.  The  model 
also  incorporates  the  influence  of  methanol  crossover  from 
anode  to  cathode  based  on  a  combination  of  diffusion,  electro- 
osmotic  drag,  and  pressure.  At  fixed  current  densities,  higher 
temperatures  of  operation  increase  the  cell  potential  and 
power  density  values.  On  scale-up,  the  cell  performance,  at 
the  moment,  tends  to  be  slightly  inferior  and  this  is  being 
explored  further.  In  conjunction  with  this,  research  is  now 
under  consideration  to  reduce  the  extent  of  methanol  cross¬ 
over  and  improve  the  performance  of  the  anode  catalysts. 
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Appendix  A.  Values  of  parameters  used  in  the  model 


e 
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cat 


an 


“an 
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cat 


(  Vo)  cat 
£o*cat 

( Vo ) 

P'0, an 
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P 

1  cat 
^m,an 
^m.cat 

^h2o 

K 

k 

C0,  MeOH 

X 

T 

R 

^MeOH 

*MeOH 


0.3 

RT IF  V / factor  of  e 
46(77353)  mV/factor  of  e 

1.5  X  10”3  cm 

5  X  10"3  cm 

6 
4 

96  488  C/equiv. 

7.14  A/cm3  @  298  K 
0.355  V 

6.25  A/cm3  @  333  K 
0.265  V 

1  atm 

2  atm 
0.17 
0.05 

2.5  H20/H+ 

6186 exp(  -7100/T)  cm2s_1  atm"1 
5.926  X  108  exp(  -9756/T)  cm  s"1 
0.0206  cm 
1  X  10"3  mol  cm"3 
1.51  V/ (mol  cm_2s“1) 

353  K 

8.314  J  mol" 1  K"1 
2.48  X  10"2MeOH/H+ 

0.0184 


Appendix  B.  Model  of  the  diffusion  gas  layer 


The  model  for  the  diffusion  mass  transfer  is  described  here 
on  the  basis  that  all  components  are  gaseous  and  well  mixed. 
The  purpose  of  the  model  is  to  determine  the  variation  in 
concentration  of  species  in  the  porous  carbon  backing  layers 
and  thus  the  concentration  of  the  reactive  species  at  the  edge 
of  the  porous  catalyst  layers.  Diffusion  in  the  porous  backing 
layer  of  the  anode  can  be  described  using  the  Stefan-Maxwell 
equations  for  a  multicomponent  gas  mixture 


n  RJ 

v*,=  I  =  1 2’K,n  (Al) 

N  is  the  molar  flux  of  the  appropriate  subscript  species,  x  is 
the  mole  fraction  of  the  species,  and  Z>eff  is  the  effective 
diffusion  coefficient  of  the  i-j  species  pair.  Using  the  sub¬ 
scripts  1,  2  and  3  for  methanol,  water  and  carbon  dioxide, 
respectively,  and  defining 


RT 

PD% 


(A2) 


dxx 

—=ku(xxN2-x2Nx)  +  kx3(xxN3-x3Nx) 
dz 

(A3) 

dx2  , 

^21  (*2^1  —xlN2)  +  k23(x2N3  —  x3N2) 

dz 

(A4) 

=  k3X  (x3Nx  - xxN3)  +  k32(x3N2~x2N3) 
dz 

(A5) 

Noting  that  kab  —  kba,  and  that 

Xx  +*2+*3=  1 

(A6) 

and  noting  the  stoichiometric  equation  for  the  methanol  oxi- 

dation  reaction 

CH3OH  +  H20  C02  +  6H  +  +  6e  “ 

(A7) 

so  that  (assuming  no  net  water  transport 
membrane) 

across  the 

11 

a? 

i 

11 

a? 

11 

(A8) 

we  can  rewrite  (A3),  (A4),  (A5)  as 

1  dxx 

Tj~T~ +  (^13  ^12)^2  ^13 

N  dz 

(A9) 

1  dx:2 

7r~T~  =  (^23  “”^12)^1  +  ^\2X2~^23 

N  dz 

(A10) 

The  solution  to  these  equations  is  obtained 
writing 

as  follows 

Nkl2  =  a  Nkl3  =  b  Nk23  =  c 

we  obtain 

dxx 

—=ax  1  +  ( b-a)x2-b 
dz 

(All) 

,  =(c  a)xx+ax-,  —  c 
dz 

(A12) 

This  is  a  system  of  first  order  ODEs:  we  can  rearrange  as 


d2JC]  dxx 

— j  ~  2a~ — b  ( ab  —  bc-\-ca)xx  —ab  =  0 


dr 


dz 


(A13) 


a  second  order  differential,  which  can  be  solved  to  give  ( when 
(a  —  b)(a  —  c)  is  positive) 


—  „(«  +  /3)z  _j_  „  „{a  -  p)z 


xx~cxe 


+  c2e 


+  y 


(A14) 


1 


*2  = 


(b-a) 


{cx^a+^-c2^eia~^  +  ay+b)  (A15) 


*3 


—  C\€^a+^z 


<ai6> 


we  can  write 


where 
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P  =  -/(a-b)(a-c) 
ab 


y— 

(ab  —  bc  +  ca) 

f  1  \ 

Xio-y-~((a~-b)x20  +  ay+b) 

\  P  ) 

(  1  \ 

xio-y+-((a-b)x2o+ay+b) 

\  p  / 

As  a  first-order  approximation  it  is  possible  to  generate 
linearised  versions  of  the  equations  which  enable  order  of 
magnitude  estimates  of  the  impact  of  diffusion  on  the  varia¬ 
tion  of  concentration  in  the  porous  backing  layers.  This  pair 
of  equations  can  be  expanded  in  a  Taylor  series  around  z  —  0 
by  differentiating  wrt  z  and  setting 

*i(0)=*i0 


Table  Al 


Values  of  g for  binary  systems  in  a  DMFC 

h2o 

o2 

CH3OH 

02 

ch3oh 

co2 

n2 

4.243  X  10_1 
1.991  X10“7 
2.731  X10"7 
4.439  X10"7 

2.759  X 10“  6 

4.367X10“* 

2.059X10"6 

2.768X10-6 

Using  the  Slattery-Bird  formulae  to  calculate  values  of  PDab> 
we  have 


PDab~  Sab6? 


2,334 


for  systems  containing  water,  and 


*2(0)  “ *20 

thus,  to  the  first  term  in  z, 


*1  U)  -*10  +  ^z(^12*10+(^13-^12)*20-^13)  (A17) 

X2(z)=*20  +  Az((^23”^12)*10  +  ^12*20-^23)  (A18) 

and,  applying  Eq.  (A6), 

*3 ( Z)  =*30-AZ(^23*10  +  ^13*20-  (kl3+k23))  (A19) 

This  approximation  is  valid  for  small  products,  Nz(kab).  As 
a  further  approximation,  in  the  special  case  where 
k\2  ~  k\i =  ^23 =  k  we  can  write 

*i  (z)  -x10  +  Az/:(x10-  1) 


x2(z)  =X20  +  Nzk(x 20-  1) 

The  value  of  Nz{kab) 

To  evaluate  the  product  Nz(kab)  we  can  write 


RT 

PD<xb 


1 


ld  =  0.25  mm 


The  PDab  product  can  be  estimated  using  the  Slattery  and 
Bird  correlation  [  19] 

1.832 

( PcaPcb)m(TcJcb)sn 2 


PDab = 0.0002745 


/  T  ; 

Wrcrf; 


X 


€ 


where  e  represents  a  correction  factor  accounting  for  the 
porosity  and  tortuosity  of  the  porous  medium:  e  is  commonly 
found  to  be  between  0.25  and  0.35. 

In  the  case  of  binary  systems  containing  water,  the  expres¬ 
sion  is 


PDab  =  0.0005 148 


(pCaPb)l/3(TCaTl) 


5/12 


X 


€ 


PDab  =  gab*T 


1.832 


otherwise.  Values  for  gah  are  given  in  Table  Al . 

Using  a  typical  value  for  e  of  0.3,  and  a  temperature  of 
80°C,  the  lowest  PDab  value  is  0.05281  for  water/ methanol. 
With  a  diffusing  region  0.3  mm  thick,  the  Nkl  product  for 
this  system  is 

RTIld 

- i  =  2.88Xl0~3/ 

6  FPD 

where  I  is  the  current  density  in  A  cm-2.  This  implies  that 
for  current  densities  less  than  3  A  cm  “ 2,  the  partial  pressures 
at  the  catalyst  layer  will  equal  the  partial  pressures  in  the  flow 
channels  with  an  error  of  less  than  1%. 

A  similar  argument  applies  to  the  cathode  diffusion  region. 


Notation 


c 

C\ 

c2 


Ac 

F 

P 

k 

K 


permeation  rate 

effective  diffusion  coefficient  of  i  in  the 

membrane 

concentration 

concentration  at  feed/ membrane  interface 
concentration  at  permeate/membrane 
interface 

C2-CJ 

Faraday  constant 
membrane  thickness 
pressure 

rate  of  permeate  removal 

constant  related  to  effective  hydraulic 

permeability 
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X, 

R 

P 

mff 

ij 

Nt 

z 


F 

T 

4i 

superscript  c 
l 


mole  fraction  of  species  i 

gas  constant,  8.3143  J/ (mol)  (K) 

pressure,  atm 

effective  binary  diffusion  coefficient  for 
species  i,j,  cm2/s 
superficial  gas  flux  of  species  i, 
mol/ (cm2)  (s) 

direction  normal  to  surface  of  electrode 
mole  fraction  of  species  i  in  the  gas  flow 
channel 

Faraday’s  constant,  96,488  C/eq 
temperature 

thickness  of  gas  diffusion  layer,  cm 

critical  value 

current  density,  A  cm-2 
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